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Using gauge/gravity duality, we study the thermalization of strongly coupled A/" = 4 supersymmet- 
ric Yang-Mills plasma. We analyze the expectation value of the stress tensor and scalar correlation 
functions and the applicability of the fluctuation dissipation theorem. Via gauge/gravity duality, 
this maps into studying the equilibration of a black hole geometry and its Hawking radiation. 



Introduction. — Gauge/gravity duality [Tj, or hologra- 
phy, is a powerful tool for studying real-time dynamics in 
strongly coupled quantum field theories. Through holog- 
raphy the creation and thermalization of non-abelian 
plasma maps onto the process of gravitational collapse 
and black hole thermalization, which can be studied nu- 
merically. The equilibration of holographic plasma can 
yield valuable insight into the thermalization of quark- 
gluon plasma produced in heavy ion collisions at the Rel- 
ativistic Heavy Ion Collider and the Large Hadron Col- 
lider. The simplest theory to study with a gravitational 
dual is A/" = 4 supersymmetric Yang-Mills theory (SYM) 
at large and 't Hooft coupling A. 

A necessary condition for a system to be thermalized 
is that it has a well defined temperature T. In a non- 
equilibrium state the SYM proper energy density £{x) 
yields one definition of an effective temperature 
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Another definition of an effective temperature comes 
from two point functions g{x\x'). Let gsymix\x') = 
^{{d{x),6{x')}) and g,nti{x\x') - -z{[6{x),d{x')]) be 
symmetrized and antisymmetrized correlation functions 
of a local bosonic operator 0{x). In a fully equilibrated 
system the Fourier transforms of these two correlators are 
related by the Fluctuation Dissipation Theorem (FDT) 
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where lo = . In a non-equilibrium system we take the 
Wigner transform 



{x,q)^ / dAxg(x|x')e-^«•^^ 
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with Ax = X — x' and x — ^(a; -|- x'), and define an 
effective temperature based on the FDT 
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In non-equilibrium states Ti{x) ^ T2{x,q) and the FDT 
is not satisfied. However, as the system thermalizes in the 
vicinity of the point x, we must have T2{x, q), Ti{x) — >■ T 



and equivalently, the applicability of the FDT. When 
this condition is satisfied at momentum q, by for exam- 
ple current-current correlation functions, the plasma will 
emit a thermal flux of photons. Thus the applicability of 
the FDT provides one of the best measures of thermal- 
ization. 

We are interested in exploring the physics of thermal- 
ization in the simplest setting which allows for complete 
theoretical control. This leads us to focus on the dynam- 
ics of homogeneous, but anisotropic, states. To create 
states such states, it is natural to consider the response 
of an initially equilibrium state in SYM to a temporally 
localized time-dependent change in the spatial geometry 
[2]. Following [2] we limit attention to the metric 

ds^ = ^dv^ + dxl + e-2b(") dxf^ , (5) 

where xj_ = {xi,X2} and v is time. For fe'(w) 7^ 0, the 
changing geometry does work on the quantum system 
and excites the state. 

In the dual gravitational description the metric ([s]) 
acts as a source for Einstein's equations, creating a non- 
equilibrium black hole via the emission of gravitational 
radiation from the time-dependent boundary geometry 
of asymptotically AdSs spacetime . There are two dis- 
tinct processes by which the black hole thermalizes. The 
first is classical and involves the relaxation of the black 
hole geometry to local equilibrium. The second process is 
quantum mechanical and involves the relaxation of quan- 
tum fluctuations near the black hole's event horizon to 
equilibrium. These quantum fluctuations — Hawking ra- 
diation — are encoded in bulk correlation functions. The 
near-boundary behavior of the gravitational field encodes 
the SYM stress T'^'^ p] and hence the effective tempera- 
ture Ti . The near-boundary behavior of bulk correlation 
functions encodes SYM correlation functions and hence 
the temperature T2. Therefore, to get Ti and T2 one 
must solve Einstein's equations and then compute bulk 
correlation functions on top of the time-dependent geom- 
etry. T2 can be regarded as the effective temperature of 
the black hole's Hawking radiation at the boundary. 

Our numerical scheme for solving Einstein's equations 
is outlined in [T, T . Diffeomorphism and translation in- 
variance allows one to chose the 5d metric ansatz 
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where A, B, and E are all functions of the AdS radial 
coordinate r and time v only. The coordinates v and r are 
generalized Eddington-Finkelstein coordinates. Infalling 
radial null geodesies have constant values of v (as well 
as x± and Outgoing radial null geodesies satisfy 

dr/dv — ^A. The boundary of the geometry is located 
at r = oo. Demanding that the boundary metric is that 
of ^ equates to imposing the boundary conditions [2 
limr->oo B{v, r) = b{v), lim^^oo ^{v, r)/r = 1. 

We choose b{v) = c/(V2^)e-^'/(2^'). As an IR regu- 
lator, we choose initial data corresponding to an equilib- 
rium black brane with temperature Tj in the infinite past 
when b{v) = 0. Once Einstein's equations are solved, 
the SYM stress-energy tensor can be extracted from the 
near-boundary behavior of the 5c? gravitational field [3]. 
We adjust the width and amplitude of b{v) such that the 
temperature of the SYM plasma in the infinite future 
(again when b{v) — 0) is Tf = SO^/^Ti. Our choice of 
amplitude and width of b{v) are c — 2.11, a — 1/ttT{. 
With these parameters the SYM energy density changes 
by a factor of 50 over a time l/irTf. 

For simplicity, we limit our attention to bulk dilaton 
correlators. The near-boundary asymtototics of the dila- 
ton correlators encodes corresponding SYM Lagragian 
correlators. With the SYM correlators known, we com- 
pute the effective temperatures Ti and T2 and analyze 
how the system thermalizes as a function of time and 
momentum. 

Bulk correlators. — Our calculation of bulk correlators 
largely follows that of where string correlators were 
considered. Because the bulk geometry is translationally 
invariant, it is convenient to introduce a spatial Fourier 
transform and compute correlators as a function of v, 
momentum q and r. The symmeterized and antisymme- 
terized dilaton correlators (both denoted generically with 
the shorthand notation G(l|2) = G{vi,ri\v2,r2) with the 
momentum dependence suppressed) satisfy the equation 
of motion 



Di^cm = -Df.^cm = 0, 
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where £>^„) = gMND^n)^^n) ^i^h D^^^ the covariant 
derivative operator with respect to coordinate label n 
(no sum on n implied). In Eddington-Finkelstein co- 
ordinates Z3? , is first order in time derivatives, so ini- 
tial data require to solve ([7]) consists of G°(l|2) = 
G(l|2)|t,j^=^o^^2=i,o. Given initial data, the solution to 
reads [I^ 



G(l|2) 
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where the retarded Greens function satisfies 



^(i)Gr(1|2) 
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and g = det gMN- At the boundary the retarded Greens 
function satisfies the boundary condition of no non- 
normalizable modes turned on. 

The initial data G°, which is specified in the past when 
the geometry was a static black brane, cannot be ar- 
bitrary. It must contain coincident point singularities, 
which are universal in form. At fixed spatial momentum 
q, the singular part of the initial data reads 

G°_ ^ , log |(«i-«2)(ri-r2)|, (10a) 
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sgn(tii-U2) - sgn(ri-r2)] , 
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where K — 2^1 with k| the bd gravitational constant. In 
the limit w° — >■ —00, only the singular part of the initial 
data contributes to the evolution of the bulk correlators 
[5]. The retarded Greens functions in ([8| propagate ini- 
tial data up to the boundary, where it reflects and gets 
absorbed by the horizon at r = r^,. This generically hap- 
pens in a time ~ 2/7rTi. However, initial data arbitrarily 
close to the event horizon takes much longer to reach the 
boundary. In the limit w° ^ — l/TrTj, only initial data 
exponentially close to the horizon contributes to the evo- 
lution of the correlators near w = 0, when the geometry 
is changing. Moreover, due to the gravitational redshift 
of the black hole, only initial data with divergent radial 
momentum will give rise to a finite wavelength excitation 
in the bulk near w = 0. 

As discussed in ^ , the above argument suggest a strat- 
egy for computing bulk correlation functions. First, in 
the limit ti° <C — I/ttTi all initial data except that expo- 
nentially close to the horizon can be neglected. Initial 
data exponentially close to the horizon can be evolved 
until it reaches the stretched horizon at r* = rh + e, 
where it determines an effective horizon correlator. The 
horizon correlator acts as a source of radiation on the 
stretched horizon, which subsequently propagates up to 
the boundary and determines the SYM correlators. 

Let Qyi be a retarded Greens function in the region 
rh < r < r^, which satisfies the boundary condition 
t?R(l|2) = at r„ = r*. The composition law for Greens 
functions, which says how information at a point 3 inside 
the stretched horizon is propagated to a point 1 outside 
the stretched horizon, reads 



(11) 



Gfl,(l|3) = - j GR(l|2)i?[2)gR(2|3). 

Substituting (111 into (|8|, we find 

G(l|2) = j Gr(1|1')Gr(2|2')G''(1'|2') 

where the horizon corerlator is 

Gni|2)=i^[i)i^[2)5(l|2)L^^,^^,^, 
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and G is given by ([8]) with the replacements Gr 

and with all initial data except that exponentially close 

horizon neglected. 

In the limit e — >■ 0, C/ can be computed using geometric 
optics. This is because the radial momentum of G must 
be divergent when e — > 0. At leading order the geomet- 
ric optics approximation leads to the equation of motion 



convolution integrals in (12 1 in frequency space, where 



(„)C+(„) 



0, where d- = dr is 



the directional derivative on infalling null radial geodesies 
and d+ = dy + ^Adr is the directional derivative along 
outgoing null radial geodesies. The general solution to 
the above equation which satisfies Dirichlet boundary 
conditions at the stretched horizon reads 
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/(Z(1),Z(2))-/(Z(U),Z(2)) 



-/(Z(1),Z(2,)) + /(Z(U),Z(2,)) , (14) 



where the function / is determined by initial data, Z 
is constant on outgoing null radial geodesies, and the * 
subscript implies the corresponding radial coordinate is 
evaluated at r = . 

Near the horizon Z = {r — rii{v))exp ~ dv' k{v') 
where k{v) = ^drA{r,v)\r=rh- By taking the coincident 



point limits v° 



V2 — ^ v° and matching the solution 



(14) onto the initial data (10), we find 
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/sym(2:i,X2) = - — log |xi - X2I, (15a) 
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Combing Eqs. (13), (14) and (15) we secure 
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X < k{v) n{v')aic\Y' I dv k{v) \ , 



gt,,{v,r\v'y) = 2K 



5'{v-v'). (16b) 



These formulas are valid at leading order in e. It is possi- 
ble to resum all of the momentum-dependent corrections 
to the horizon correlators, which are in powers of eq^ . 

The metric ^ is invariant under the residual diffeo- 
morphism r — > r -I- S^{v) where S^{v) is arbitrary. After 
solving Einstein's equations and finding the horizon of 
the geometry, we exploit this residual diffeomorphism 
invariance to set the horizon position — 1. Choos- 
ing e = 0.015, we then numerically solve ^ for the re- 
tarded Greens function Gr, compute the horizon cor- 



relators via (16), and compute the bulk correlators by 



numerically evaluating the integrals in (12). The sym- 
metrized horizon correlator diverges like l/(w — w')^ a-t 
coincident points. Therefore, it is convenient to do the 



the coincident point singularity can be analyzed analyt- 
ically. To increase the rate of convergence of the fre- 
quency integrals, we choose to compute the smeared cor- 
relators G(w|w') = / dvidv2huj{v — vi)hw{v' — V2)G{vi\v2) 
where h^{v) is a Gaussian of width w — 1/3ttT{. With 
the smeared bulk correlators known, the corresponding 
smeared SYM correlators are given by [5] 



g{v\v') = (16/if2) lim 



r^r'^G{v, r\v' , r 



(17) 



We extract the effective temperature T2 from the 
smeared SYM correlators via Eqs. ([s]) and ([4|. We have 
checked the e dependence of our numerical results by 
changing e by a factor of two. This yields a 3% change 
in our results for the smeared correlators, indicating con- 
vergence. We comment more on the numerical details of 
our calculation in a coming paper. 

Discussion. — Fig. [l|a) shows a plot of the SYM stress 
T^\, = diag(— £',Pj_,7^j^,7^||) as a function of time v. In 
the distant past the stress is static and corresponds to a 
low temperature equilibrium plasma of temperature Tj. 
When the boundary geometry ^ starts to change, work 
is done on the system, the energy density generally grows 
and the pressures oscillate. The energy density becomes 
constant when the boundary geometry becomes static, 
and is within 2% of its final equilibrium value at time 
V ~ 1.5/ttT[. The transverse and longitudinal pressures 
and 1 equilibrate near time v = S/nTf. 

Fig. [ijb) shows the effective temperatures Ti{v) and 
T2{v, q) for time-like momentum (a; = TrTf , q± = q\\ =0), 
and light-like momentum (a; — 7rT{,qj_ = = 7rTf/-\/2), 
where q± and are the spatial momentum in the trans- 
verse and longitudinal directions. In the distant past, 
the system is in equilibrium and the effective tempera- 
tures are static and identical. When the boundary ge- 
ometry becomes dynamic, the effective temperatures all 
differ from each other. Moreover, T2 is singular during 
the non-equilibrium evolution. Clearly it makes no sense 
to interpret the effective temperatures as true tempera- 
tures during the dynamic stage of the evolution since they 
do not agree with each other and can become singular. 
However, as time progresses the effective temperatures 
all asymptote to the same constant value Tf . 

Examining Fig. [l|b), we see that Ti approaches Tf 
before T2. Furthermore, the time- like mode of T2 ap- 
proaches Tf before the light-like mode. This behavior is 
easy to understand from the dual gravitational physics. 
In the dual gravitational description, the equilibration 
of the SYM stress reflects the equilibration of the near- 
boundary geometry. The near-horizon geometry equili- 
brates as non-equilibrium gravitational modes fall from 
the boundary through the horizon. The dilaton hori- 
zon correlators equilibrate at the same rate as the scalar 
components of the near-horizon metric (which in our 
setup equilibrate faster than the anisotropic components 
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oo. In this limit one can again use geometric 



optics to compute the retarded Greens functions in (12 1 
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FIG. 1: (a) The SYM stress tensor as a function of time, (b) 
The effective temperatures Ti and T2 as a function of time. 
T2 is plotted for time-like momentum uj = nTf and q — and 
light-like momentum lj = yrTf and q± = = 7vTf/\^. Ti 
vanishes when £ = and T2 vanishes when 2gsym = iipanti 
and diverges when ganti = 0. 



of the metric). When this happens the horizon corre- 
lators satisfy a horizon fluctuation-dissipation relation, 
(i.e. Eq. ([2I with g ^ Q^) and the flux of Hawking 
radiation from the horizon is thermal. However, an ob- 
server near the boundary will still see a non-equilibrium 
flux of Hawking radiation. For generic mode q, it takes 
an additional time Av > 2/7rTf for the thermal flux of 
Hawking radiation to propagate up from the stretched 
horizon to the boundary, and hence for mode q of the 
SYM correlators to thermalize. 

The dual gravitational physics can also tell us which 
modes thermalize first. The equilibrated modes which 
first arrive at the boundary from the stretched hori- 
zon propagate along null radial geodesies with a transit 
time of order Aw ~ 2/TrTi. These geodesies are relevant 
for the geometric optics approximation of the retarded 
Greens functions in ( 12 ) in the limit a; — > 00 with q 
fixed. Hence the modes that thermalize the fastest have 
high frequency and small momentum. Similar conclu- 



The relevant geodesies travel very far in the spatial di- 
rections and at asymptotically late times approach the 
boundary as slow as r = (7rTf)^z;. However, for oj = \q\ 

the geometric optics approximation breaks down when 

1/3 

r - TrTf (|q|/7rTf)"'^ We therefore see that the thermal- 
ization time for large momentum light-like modes can 
be as large as Av - (l/TrTf) (|q|/7rTf)^/^ In the dual 
field theory the origin of the \q\^^^ upper bound on the 
thermalization time lies in the fact that nearly on-shell 
excitations of momentum q can be created during the 
far-from-equilibrium evolution of the system. These ex- 
citations propagate at nearly the speed of light and have 



a lifetime 



11/3 



In Fig. [ijb) only glimpses of the delayed thermaliza- 
tion of high momentum light-like modes can be seen: the 
light-light mode of T2 approaches T{ slightly after the 
time- like mode of T2 approaches Tf . As the above argu- 
ment suggests, in our numerical simulations we have seen 
this effect becomes more dramatic as \q\ is increased in 
magnitude. However, the numerics become challenging 
in the limit u ^ \q\ ^ nTf. We further address this 
interesting regime semi-analytically in a coming paper. 
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sions were also reached in [B]. 

The slowest modes to thermalize are light-like modes 
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